Analytical soft-core potentials for macromolecular fluids and mixtures 
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An analytical description of polymer melts and their mixtures as liquids of interacting soft colloidal 
particles is obtained from liquid-state theory. The derived center-of-mass pair correlation functions 
with no adjustable parameters reproduce those computed from united atom molecular dynamics 
simulations. The coarse-grained description correctly bridges micro- and mesoscopic fluid properties. 
Molecular dynamics simulations of soft colloidal particles interacting through the calculated effective 
pair potentials are consistent with data from microscale simulations and analytical formulas. 



The formulation of an accurate mesoscopic description of macromolecular fluids has been a longstanding goal in 
polymer physics. Experimentally-relevant polymer dynamics span a wide range of timescales, for which large-scale, 
long-time properties still depend strongly on the local molecular structure m. Pertinent information on structure 
and dynamics of polymer liquids has been gained from united atom (UA) molecular dynamics (MD) simulations. 
However, the MD computational time increases as the squared number of interacting units, and the latter has to be 
large to approximate the thermodynamic limit, rendering an all-atom simulation of long-time polymer dynamics a 
prohibitive task. One strategy devised to overcome this problem is to renormalize the liquid structure and dynamics 
using effective- unit coarse-grained descriptions Specifically, polymers can be described mesoscopically as soft 
interpenetrating spheres having the overall size of the polymer, i.e., the radius of gyration Rg. However, to correctly 
perform the renormalization procedure, a theoretical framework that bridges properly different lengthscales of interest 
is needed. Phenomenological mesoscopic potentials were implemented by Dautenhahn and Hall, and later on by Murat 
and Kremers, to describe polymer melts and blends ^ JjJ. Hansen and coworkers have recently developed a rigorous 
numerical description of polymer solutions as liquids of soft interacting colloidal particles . 

In this Letter we start from first-principles liquid-state theory and derive an analytical form of center-of-mass 
(c.o.m.) pair correlation functions, from which the effective pair soft-core potential acting between molecules in 
polymer liquids (melts) and their mixtures (blends) is obtained. The c.o.m. pair correlation functions reproduce 
mesoscale liquid structures obtained from UA-MD simulations 0,0.13 without adjustable parameters. Test systems 
are polymer melts with different architecture, local semiflexibility, and degree of polymerization (Table I), as well as 
their mixtures (Table II) . Finally, the mesoscopic potential derived by an inversion procedure is used in MD simulations 
of soft colloidal particles, which reproduce the liquid structure at the level of c.o.m. pair correlation functions. 

The renormalized pair interaction potential is a function of the c.o.m. total pair correlation function h{r). In 
reciprocal space, h{k) = [cj™'^(fc)/[x)"""(A;)]^ft,"""(fc), after a procedure devised by Krakoviack, Hansen, and Louis 8]. 
Here, uj"^'~' (k) is the intrachain monomer distribution about the c.o.m. Also, h™'™{k) and uj^"-™{k) are the monomer- 
monomer intermolecular total pair and intramolecular correlation functions, respectively. We assume a Gaussian 
description for intramolecular form factors: w"""(A:) = N/{1 + k^R'^/2), where N is the number of "monomeric" 

units of length a and Rg = ^N/&a is the overall chain size, and a;™'^(fc) — Ne~^'^^^!il^ For h™'™{k), we use the 
thread-limit PRISM description in which the chain is treated as an infinite thread of vanishing thickness. Since 
this model correctly reproduces the appearance of correlation hole effects on the scale oi Rg, a distinguishing feature 
of macromolecular fiuids, liquid properties at the mesoscale of interest should be described well. Now, /i™™(fc) — 
Att^'p [Cp/(1 + fc^^p) - Cc/(1 + ^^C?)] I where = Rg/V^ is the correlation hole lengthscale. The density fluctuation 
length is deflned as = +Cp'^ with = Rg/{2np*i^), and p*^ — PchR^, where pch is the chain number density. 
In real space, we find 
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where erfc(a;) is the complementary error function. This expression satisfies the condition h{r) > — 1, spanning melt 
to dilute solution densities up to p*^ > 0.03. Eq. Q correctly recovers the isothermal compressibility kt E3 related 
to the k ^ hmit of the mesoscale static structure factor as S{0) = 1-1- pchh{0) = (Cp/Cc)^ = PchkBTniT- 

Our total correlation function h(r) effectively maps the polymeric liquid onto a fluid of soft interpenetrating colloids 
of radius Rg and density pch- Translated into the colloidal-particle framework, Eq. is a function of the reduced 

variables f = r/Rg, Cp = Cp/Rg = (27rp*;,)"\ and = ^p/Rg = [%/2(l + 7r%/2p*^)] "\ Then, Eq. ^ reduces to 



2 



h{r) = /i(r, ^p), where we define 
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We tested Eq. ^ against UA-MD simulations for polymer melts of polyethylene (PE) for increasing degree of 
polymerization {N — 44, 66, 96) ^ as well as different local structure and flexibility, including syndiotactic (sPP), 
isotactic (iPP) and head-to-head (hhPP) polypropylenes |EI3- Input parameters to the theory were site number 
density p, temperature T, and Rg from Table I. Our analytical expression agrees well with h(r) from simulations within 
statistical error (Figs. 1 and 2). Specifically, for PE melts it reproduces the tendency for chains to interpenetrate 
more efficiently with increasing length, an effect due to the fractal dimension of polymer chains. For slightly branched 
polymer melts, the number of interpolymer contacts becomes larger with increasing polymer stiffness (increasing Rg 
at constant N). Good agreement between theory and computer simulations is found also in reciprocal space, where 
the analytical Fourier transform of Eq. ^ reproduces simulation data in the entire range for k > R^^ (Fig. 3). 

A simplified form of Eq. can be derived in the limit of long polymer chains (fp — > 0) as is implicitly assumed by 
the use of Gaussian form factors, for which we find 



where h{f,£^p)/$,p is a universal function of the reduced distance f up to the second-order correction in ^p. Eq. © 
is a good approximation for Eq. (Q): results obtained with the two expressions cannot be distinguished in Fig. 1. 
Both equations recover the correct trend of greater chain interpenetration with increasing Rg and/or liquid density 
asMO)(x (p:J-i. 

From Eq. {Q), the effective potential is calculated self-consistently using an inversion procedure based on the HNC 
approximation, (3v{r) = h{r) — \n[h{r) -I- 1] — c(r), where (3 = l/ksT. The direct correlation function c(r) ll] is 
obtained from the Fourier transform of the Ornstein-Zernike equation c(k) — h{k)/ [1 + pchh{k)]. While the HNC 
closure is known to work well for dilute colloidal systems, its applicability is questionable in dense systems where 
many-body interactions should be important pT| . However, the pair HNC interaction potential works well for the 
mesoscopically-renormalized polymer melts investigated here. We performed constant temperature MD simulations 
of a liquid of soft colloidal particles of radius Rg interacting through the derived pair potential, for chains of = 44, 
66, and 96 units. Corresponding effective potentials Pv{r) as functions of the normalized distance r/Rg differ slightly 
for the three cases. The potential for A = 44 is shown in the inset of Fig. 1. Mesoscale simulations recover the liquid 
structure obtained in UA-MD simulations in both real (Fig. 1) and reciprocal space (Fig. 3). 

We extended the formalism just presented to treat binary polymer mixtures. The homopolymer species are 
A and B, characterized by radii of gyration RgA and RgB, with number of chain units Na and Nb, and unit 
lengths a A = ^jQ/NARgA and as — \/Q/NBRgB- The volume fraction of component A is given by (j) and 
7 = aB/o'A represents the mismatch in local chain semiflexibility. For a generic pair {a,P £ (A, i?)}, we have 
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hap{k) = Ci;™'^(fc)w™'-^(fc)/(aj™™(A:)ci;™™(fc)) h'^™{k). Chains are assumed to obey Gaussian intramolecular distri- 
butions with uj''^'^{k) and w""(A:) defined as in the melt case. The /i™"(/c) follow the PRISM-blend thread model 
described by Tang and Schweizer jl^ . and is extended here to include asymmetries in local chemical structure and 
fiexibility {a a < cb). The advantage of this approach is that no closure approximations need to be invoked, while the 
miscibility parameter x enters as an input to the theory [l^ . This allows us to describe the renormalized structure 
and potential for systems following cither upper or lower critical solution temperature phase diagrams within one 
theoretical framework. 

On applying inverse Fourier transforms, the blend haf}{k) in real space are given by 

hABir) = -lUr)+I'ABir) , (4) 

with 
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Here for {A G (p, 0)} identifies the lengthscale for density or concentration fluctuation correlations, 
concentration fluctuation length £^ — a ab I \/24(:Vs — y)<^(1 — 4)) diverges at the spinodal temperature where 



length = ^TAB/y24(^ 
The average segment length is ct^^ = (f)c 

with being the average correlation hole length. Also 
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The self terms in the density fluctuation contributions I^q, (r) 

are formally identical to the total pair correlation function in the melt. For a totally symmetric blend {Na = Nb with 
CTA = (^b), Eqs. (|4I5|) correctly recover the melt equation in the athermal limit where enthalpic effects are negligible. 

We tested Eqs. (|4I5|) against UA-MD simulations of blends with polymer components having different degree 
of polymerization N, local semiflexibility /branching cr, and volume fraction 0. By convention, the i?-component 
is assumed to be the stiffest. In Fig. 4, we report only two representative examples (see Table II). For a thermal 
(^-symmetric mixture of linear (PE) and slightly branched (hhPP) molecules, the theory agrees well with simulations 
(x = 0.0016 Ref. 13] and Xs = 0.0211). The stiffest component (PE) shows a higher number of intermolecular 
contacts than the flexible one (hhPP). A weak self-clustering of the most flexible species is observed at a distance 
comparable to the polymer size, while the number of AB contacts is still large at this temperature, an indication that 
the system is far from its phase transition. Similar effects are shown for the second system in Fig. 4, which is a thermal 
0-asymmetric blend of slightly branched (iPP) and linear (PE) molecules (x — 0.005 Ref. and Xs — 0.0281). In 
asymmetric mixtures, the minority species (iPP) presents an enhanced clustering, in agreement with experimental 
data and simulations 0, Q ■ 

The blend pair correlation functions can be mapped onto a system of interacting soft colloids of dimension RgA and 
RgB , with chain volume fraction given by = (f>ch and mismatch in chain size given by 7 = Rgs / RgA ■ The analytical 
equations are considerably simplified for symmetric mixtures where 7 = 1, and we confine our presentation to this 
case Q. For each pair interaction, Eq. Q reduces to Eq. I^ir, ix) = h{f, ^x) where ^x for {A € (p, 0)} and — 
£.<p/Rg = l/\/2(l — x/Xs)- By combining Eqs. 14I5|I into density and concentration fluctuation contributions following 
Bhatia-Thornton's formalism we recover known expressions for a mixture of symmetric colloidal particles (e.g., 
liquid alloys). For example, in the k limit, the correlation in number density S'^p{0) — (Cp/Cc)^ correctly recovers 
kt, while the concentration fluctuation correlation reduces to the known formula S'^'^{0) = — (f>)/[l — 2(f){l — (f>)Nx] 
upon introducing Flory's definition of Xs- Finally, the cross term 5'''*(0) vanishes in agreement with the theory of 
symmetric particle mixtures Consistency between Eqs. H4I5|) and known properties of colloidal mixtures supports 
our mapping of polymer blends onto a mixture of interacting soft colloidal particles. 

Summarizing, we report here the derivation of an analytical renormalized description of polymer melts and blends 
as fluids of mesoscopic soft colloidal particles. The related soft-core Gaussian potential explicitly bridges meso- and 
microscale properties. The derived c.o.m. pair correlation functions reproduce, in both r- and A;-space, fluid structures 
obtained from UA-MD simulations. Using the melt potential calculated from the HNC approximation, we perform MD 
simulations of liquids comprised of soft colloidal particles, from which the obtained pair correlation functions correctly 
recover the liquid structure from UA-MD simulations, further supporting our renormalized analytical description of 
polymer liquids. 
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TABLE I. Polyolefin melts. 



Polymer 


N 


T[K] 


p [sites/ A-'^] 


R, [A] 


PE 


44 


400 


0.0324 


10.50 


PE 


66 


448 


0.0329 


13.32 


PE 


96 


453 


0.0328 


16.78 


sPP 


96 


453 


0.0328 


13.93 


hhPP 


96 


453 


0.0336 


13.54 


iPP 


96 


453 


0.0328 


11.34 



TABLE II. Polyolcfin blends {T = 453 K and Na = Nb = 96). 



Blend [A/B] 





p [sites/ A^] 


RgA [A] 


7 


hhPP/PE 
iPP/PE 


0.50 
0.25 


0.0332 
0.0328 


12.32 
11.35 


1.34 
1.47 
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FIG. 1. Plot of h{r) against r/Rg for PE melts. Theory (full lines) is compared with UA-MD (filled symbols) and 
mcsoscale (open symbols) simulations for iV = 96 (squares), 66 (circles), and 44 (diamonds). Inset shows a plot of 
Pv{r) against r/Rg for the iV = 44 case. 

FIG. 2. Plot of h{r) against r/Rg for melts of different polyolefins. Theory (full lines) is compared with UA-MD 
simulations for PE (squares), sPP (circles), iPP (downward triangles). Inset shows hhPP (upward triangles). 

FIG. 3. Plot of h{k) against k > Rg^ ■ Theory (full lines) is compared with UA-MD simula- 
tions (filled symbols) and mesoscale simulations (open symbols). Left panel shows PE melts with N = 
96 (squares), 66 (circles), and 44 (diamonds). Right panel displays melts of different polyolefins: PE (circles), sPP 
(squares), iPP (downward triangles); inset shows hhPP (upward triangles). 

FIG. 4. Plot of h{r) against r for blends. Theoretical curves in athermal (solid lines) and thermal (dashed lines) 
limits are compared with AA- (circles), AB- (diamonds) and BB-teims (squares) of h{r) from UA-MD simulations 
of representative (^symmetric (upper panel) and ^- asymmetric (lower panel) blends. 



